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We obtain a lower bound on the sum of two orthogonal spin component variances in a plane. 
This gives a novel planar uncertainty relation which holds even when the Heisenberg relation is 
not useful. We investigate the asymptotic, large J limit, and derive the properties of the planar 
quantum squeezed states that saturate this uncertainty relation. These states extend the concept of 
spin squeezing to any two conjugate spin directions. We show that planar quantum squeezing can be 
achieved experimentally as the ground state of a Bose- Einstein condensate in two coupled potential 
wells with a critical attractive interaction. These states reduce interferometric phase noise at all 
phase angles simultaneously. This is useful for one-shot interferometric phase-measurements where 
the measured phase is completely unknown. Our results can also be used to derive entanglement 
criteria for multiple spins J at separated sites, with applications in quantum information. 
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I. INTRODUCTION 

Heiseiiberg's famous uncertainty relation for angu- 
lar momentum takes the well-known product form 
AJxAJy > |(iz)|/2. If (Jz) = 0, as in a singlet 
state, the Heisenberg relation becomes trivial: it only 
constrains the spin variances to positive values. Never- 
theless, there are still fundamental limits to these uncer- 
tainties, which are directly related to quantum limits on 
interferometric phase measurements [iHal) together with 
problems like macroscopic entanglement |6H9|, Bohm's 
Einstein-Podolsky-Rosen (EPR) paradox [l^l and steer- 
ing [iMi- 

Reduction of quantum noise in one spin component - 
or single-axis squeezing - is a valuable tool for enhanc- 
ing the sensitivity of interferometers and atomic clocks 
P, It has been recently implemented for ultra-cold 
atomic Bose-Einstein condensate (BEG) interferometers 
This type of quantum noise reduction reduces 
the measurement noise near some predetermined phase. 
However, if the phase is completely unknown prior to 
measurement, then it is not known which phase quadra- 
ture should be in a squeezed state. 



In this paper we demonstrate that spin operators per- 
mit another type of quantum squeezing that we call pla- 
nar quantum squeezing (PQS), which simultaneously re- 
duces the quantum noise of two orthogonal spin projec- 
tions below the standard quantum limit of J/2, while 
increasing the noise in a third dimension. This allows 
the prospect of improved phase measurements at any 
phase-angle. PQS states that reduce fluctuations every- 
where in a plane have potential utility in 'one-shot' phase 
measurements, where iterative or repeated measurement 
strategies cannot be utilized. 

Planar quantum squeezing is related to a planar un- 
certainty relation, true for any quantum state. This has 
the form of a lower bound on the planar spin variance 
sum 0,[ll, 



A'J|I > Cj . 



(1.1) 
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Here J|| is the spin projection parallel to a plane, so in 
the X ~Y plane, A^J,, ee A'^Jx + A'^Jy, and Cj is 
the minimum of the uncertainty sum for quantum states 
with fixed spin J. While values of Cj for J = 1/2 (l5| 
and J = 1 [tI were known previously, we calculate Cj 
for arbitrary spin quantum number J. We show that the 
Cj uncertainty principle has a fractional exponent be- 
havior, with Cj ~ J^^^ , and that states saturating this 
uncertainty condition have variances with the same frac- 
tional power law exponents in two orthogonal directions. 
These also have a mean spin vector in the direction of 
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Figure 1. (Color online) Three-dimensional plot of uncertain- 
ties for planar quantum squeezing. The figure corresponds 
to 100 atoms { J ~ 50) in the ground-state of a double-well 
Bose-Einstein condensate with critical attractive interaction. 
Spin variances are reduced in both axes parallel to the X — Y 
plane, with a maximum variance reduction in the mean spin 
direction X. The variance must increase perpendicular to the 
squeezing plane, along the Z axis. 



greatest variance reduction. The variance perpendicular 
to the squeezing plane is increased, with A^Jl J'^^^. 
The overall three-dimensional variance has an ellipsoidal 
shape, graphed in Fig. 1. We show that the ground state 
of a two-mode Bose-Einstein condensate (BEC) with at- 
tractive interactions gives the maximum possible PQS, 
making this an important candidate for interferometric 
measurements. 

As well as allowing improved interferometric phase 
measurements, planar uncertainty relations are useful for 
the detection of non-classical behavior in mesoscopic sys- 
tems with large total spin J. The general form of the cor- 
responding two-spin Local Uncertainty Relation (LUR) 
criterion to detect entanglement among iV sites of spin J 
particles is 



(1.2) 



where J*"*''' represents the collective spin operator, and 
the particular case of = 2 has been considered by Hof- 
mann and Takeuchi [7] as a criterion for entanglement 
between two sites. Larger N values signify multipartite 
entanglement. Experiments involving measurements in 
two spin directions [16] have employed a similar inequal- 
ity to detect genuine multipartite entanglement in four 
qubit states [l7fl . 

We extend this microscopic analysis to mesoscopic 
spins of arbitrary J, and hence show that equation (II. 2p 
is a multipartite entanglement criterion for N sites. This 



entanglement signature is valid in the mesoscopic limit of 
large J, and applies regardless of the third component, 
which may not be readily measurable. 



II. PLANAR UNCERTAINTY RELATIONS 

Our first task is to find the lower bound Cj of the 
planar uncertainty relation equation (jl.ip . We consider 
states of fixed spin dimensionality J. The most general 
pure quantum state of this type has dimension d = 2 J-l-1. 
Expressed in the Jz basis, this state is: 
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(2.1) 



Here (ft^ (m = — J, J) are real numbers indicating 
amplitude and phase respectively, and the normalization 



coefficient is 



L^ra—- 



A. Uncertainty Minimization 

To obtain a lower bound on the variance sum we min- 
imize the uncertainty over all possible expansion coeffi- 
cients, using: 



(2.2) 



Due to spherical symmetry, it is enough to consider the 
uncertainty relation in the X — Y plane. We can choose 
axes in the X — Y plane so that ( jy) = 0, with no loss 
of generality. We calculate the expectation values in the 
Z basis. We find that the squared projections are 
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(2.3) 



where we have defined J — J + 1/2. On maximizing 
the magnitude with respect to the phase variable, and 
introducing M = m -h 1/2, M± = M ± 1/2, we find that 
the mean projections satisfy: 



n ^ — ^ ^ 



-M2 



(2.4) 



Using these equations, we can numerically obtained 
the value of Cj for any spin J, by using nonlinear opti- 
mization techniques (we used the quasi-Newton method 
implemented via the Mathematica 8.0 FindMinimum 
functionfl^) to search for the minimum value of equa- 
tion (|2.2p given any possible coefficients. 

The results for selected values of J are tabulated in 
Table m and Cj/J is graphed in Fig. [Hby dots. 
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Figure 2. (Color online) Graph of planar uncertainty relation 
bound. The figure plots the lower bound to the planar uncer- 
tainty relation, normalized by J, the standard quantum limit. 
Note that C.j/J shows a decrease with increasing J. Discrete 
points are calculated numerically, while the solid line is the 
analytic approximation Cj"'. 
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Table I. Lower bound C.j of the planar uncertainty equation. 
Numerical results for equation (|l.l|l are tabulated as a func- 
tion of spin J. 



B. Asymptotic C.j values: 



We wish to obtain the asymptotic, large J limit of the 
planar uncertainty principle using an integral approxima- 
tion. The mean spin projections for any quantum state 
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m— — J 

X sin( 



{Jz) = . 



(2.5) 



For J — > oo we replace the summation over 2J-\-l discrete 
values of m by an integral over a continuous range of 
(—1,1). We define J = J + 1/2, and introduce scaled 

variables oi x — m/J and r{x) — Rm\/j, so that the 
normalization integral becomes: 



m— — J 



(2.6) 



For the mean squared spin vector, taking equation (|2.3p 
in the limit of J — cx) one obtains 
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{x)j'^ {x)dx , 



(2.7) 



where J{x) = J^/\ 

For the mean value term in equation (|2.4p . we 

also replace summations by integrals in the limit of 
J — )> cxD and define a slightly modified variable x = M / J . 
Next, using a Taylor expansion, we obtain: 



1 / N , 1 dr{x) 

r[x ± — -) — r(x) ± — - — ; — 
2J 2J dx 



8J2 dx^ ^ ^J3^ 



(2.8) 

After integrating by parts, and defining r'(x) = 
dr{x)/dx, this can be expressed as a variational calcu- 
lus problem. We minimize V [r] = A^J|| as a function of 
r{x), where 




(2.9) 
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{r'{x)) xr{x)r'{x)\ J[x)dx 
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Due to the symmetry of the integrand around x = 0, the 
minimum function r{x) must have a refiection symmetry 
with a maximum at a; = 0. In the large J limit we can 
assume a relatively narrow Gaussian solution of variance 
a <^1, where a , v is still undetermined, and: 



(2.10) 



With this choice, we can extend the integration limits to 
X = ±oo to leading order in 1/J. This means that n = 1, 
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and on expanding the integrand one finds that: 

/>CXD 

V[r]^J^ r^{x)[l- x'^]dx (2.11) 
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Next, we expand the variational integral, retaining lead- 
ing terms only. We use the results that: 



r'^{x)x^dx = 3(T^ 



(x) \/\ — x^dx — 1 



(2.12) 
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giving a corresponding asymptotic estimate for Cj of 



4a 



(2.13) 



Applying variational calculus so that dV/da = 0, and 
solving in the limit of large J, we find that a = (2J)~^^^. 
This means that: 

lim Cj w 3(2 J)2/3/8 = 0.595275 . (2.14) 

J— >oo 

Using this asymptotic form, and numerically fitting the 
tabulated results with a series in J^^^, we obtain the 
following analytic approximation to Cj: 



~ 0.595275 J^/^ - 0.1663 J^/^ + 0.0267 . 



(2.15) 



This is given in Fig. [Uby the solid curve, in good agree- 
ment with our numerical results: within 1% for J > 5. 



III. CHARACTERIZATION OF THE MINIMUM 
VARIANCE STATE 



in the Jz basis as in equation (|2.5p . As previously, we 
choose axes with squeezing in the X—Y plane and {Jy) = 
0. By selecting the phase as 0„i = 0, we get a minimum 
planar variance, and a non- vanishing mean value in (Jx)- 
This is a generic property of planar squeezed minimum 
variance states, which always have a spin vector with a 
finite mean amplitude in the plane of minimum variance. 

The value of the spin variance in each direction in the 
plane still needs to be calculated, in order to define the 
properties of the planar squeezed state. We assume a 



symmetric amplitude distribution with i?„ 



that (Jz) = 0. The mean values are then given as in 
equation (|2.4p . by: 



M=J 



(Jy) - (Jz) = . 



(3.1) 



Introducing J+ = J + 1 , r7i+ = m + 1, we find the 
squared spin projections and correlations are: 



(3.2) 
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(3.3) 
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We now wish to characterize the properties of the state 
that minimizes a planar variance sum. The main dif- 
ference can be seen if one considers that not all the 
states that minimize the Heisenberg product have pla- 
nar squeezing. For example, the state with {Jz) — —J 
is an eigenstate of the spin projection in the Z direction. 
This gives a minimum variance (zero) in one direction, 
but cannot minimize the variance sum. It is not a pla- 
nar squeezed state, since the orthogonal variances are not 
reduced below the standard quantum limit. 



A. Optimum planar squeezed state: 

We wish to analyze the detailed asymptotic properties 
of the planar quantum squeezed state that saturates the 
uncertainty principle. We consider states of fixed spin 
dimensionality J, and calculate the mean value of (J) 



(i|> = - E ' 

m— — J 



(JxJy) = (JxJz) = (JzJy) =0. 



(3.4) 



(3.5) 



We now treat the asymptotic, large J limit using an 
integral approximation as previously. For the mean value 
term, on replacing summations by integrals, in the limit 

of J — )■ 00 and defining x — M/J, r{x) = Rm\/J and 
Jx — (Jx) we obtain: 



Jx = 



2J2 



xrjxyjx) \ 
H . ^n, N — I J(x)dx, 



4J2(x) 



(3.6) 
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where J{x) = J^l - x'^. We know that a = (2 J) 
get the minimum variance. To leading order, one finds 
that: 



1 



Jx = J I r^ix) 1 



Jix)dx, (3.7) 



Next, we expand the variational integral, retaining lead- 
ing terms only, giving 



1 
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(3.8) 

Similarly, on evaluating the square of the z-projection, 
we find that: 
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(3.9) 



and hence the Heisenberg uncertainty principle in the 
Z — Y plane is obeyed, since: 



AJzAJy 



J 



>\{Jz)\/2. 



(3.15) 



On calculating the mean and variance of each spin 
component, we find that that the plane of squeezing 
always includes the mean spin direction (J). Choos- 
ing axes so that this projection is in the X direction, 
with planar squeezing in the X — Y plane, we calcu- 
late that: (Jx) ~ J - ^ ('//4)^''^ , A^Jx ~ | (2 J)^/^ , 

A'^Jy i (2J)^/^ , A^Jz - (JV2)^/^ . Asymptotically 
this is a Heisenberg limited or 'intelligent' [l^ state in 
the Z — Y plane since: 



AJyAJz « i 



(Jx) 



(3.16) 



The sum of planar variances is simple to calculate. Fol- 
lowing the techniques given previously, we find that: 

(Ji + J|.) = P(l-a)-i. (3.10) 

This means that we carry out a check on the total spin, 
which is given by the expected result of: 
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j(j + i) 



(3.11) 



However, the individual planar variances are more 
complex. Introducing scaled variables oi x = {m-\-l)/ J — 
M' J , one must use a Taylor expansion so that: 

rix ± i) - rix) ± jr\x) + i^r" {x) + 0[l^) 

. . X , . r(x) x^r(x) „ , 1 . 
= .(.)T^Kx)-^ + ^ + 0(-^). 

(3.12) 

After carrying out the integrals, we find that the squared 
projections in the plane are: 



Ji) ^ J' 
1 
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on subtracting the mean values squared, this gives the 
result that to leading order: 



A'Jx 



A2jv 



(2J) 



2/3 



(2J) 
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In summary, the important features of the optimum 
planar squeezed state are: 

• A large spin expectation parallel to the plane of 
squeezing 

• Maximum variance reduction in the direction of the 
mean spin vector 

• A smaller variance reduction in the plane of squeez- 
ing orthogonal to the mean spin 

• Saturation of the Cj uncertainty relation 

• A complementary variance increase in the third di- 
mension 

• Heisenberg limited asymptotic uncertainties per- 
pendicular to the mean spin vector 



IV. APPLICATIONS OF PLANAR QUANTUM 
SQUEEZING 

A. Generation of PQS in a BEG 



A'Jz ^ (jV2) 



We first wish to consider how to generate these PQS 
states. While there are many possible strategies, the 

(3.13) simplest is to find a physical system whose Hamilto- 
nian equals the variance sum. The ground state of such 
a Hamiltonian will minimize the variance, hence creat- 
ing a perfect planar squeezed state. This can be read- 
ily achieved in a two-mode Bose-Einstein condensate, 
which has been experimentally demonstrated to gener- 
ate spin-squeezing [3, 4]. In the limit of tight confinement 
and small numbers of atoms, this type of system can be 
treated using a simple coupled mode effective Hamilto- 
nian, where k is the inter-well tunneling rate between 
wells, and g is the intra-well s-wave interaction between 

(3.14) the atoms. The system is depicted schematically in Fig. 
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Figure 3. (Color online) (a) Schematic diagram of double-well 
Bose-Einstein condensate. Shows a schematic double poten- 
tial well, with corresponding pseudo-spin operator, (b) Spin 
variances for the ground state of a double-well BEC. Here 
A'^ = 100 atoms, or J = 50. The solid line is the total vari- 
ance, the dashed line is the Y-axis spin projection variance, 
the dot-dashed line is the X- axis variance in the mean spin 
direction. 



131 We note at that the total particle number is conserved, 
with eigenvalue N — 2J, where J is the equivalent effec- 
tive spin quantum number. 

Following standard techniques (22| - [2^ . the two- well 
BEC Hamiltonian can be written in spin language, ignor- 
ing conserved terms proportional to N or N'^. We start 
with a two-mode Hamiltonian in the standard form: 



H/h- 



K 













(4.1) 



Inter-well spin operators have already been measured 
in this environment, and are defined as: 



1 r 

2 



Jl 



{a)a-y^b^ /2, 



1 



i^b-ay^ 



(4.2) 



It is convenient to introduce a symmetry breaking vector 
Jo = {K/g,0), which causes tunneling, to give: 



H/h = -g 



II 



(4.3) 



It is clear from this form of the Hamiltonian that with 
an attractive coupling so that g < 0, the ground state will 
exactly minimize the planar variance, provided the tun- 
neling rate k is adjusted so that ^ = Jo . Since we have 

already calculated the optimum mean spin vector, we 
therefore expect that a planar squeezed state will occur 
as the ground state of the Hamiltonian for an attractive 



Jx 



1/3 



coupling {g < 0), with k/ \g\ = 

The eigenstates of this Hamiltonian can be readily cal- 
culated numerically from its matrix form. It is known to 
have a macroscopic inter-well spatial entanglement [20- 
[2^ . which is maximized at a critical attractive coupling 
value [H-lil. 

A graph of the variances against coupling is shown in 
Fig. [5]d, showing the expected reduction in both vari- 
ances at a critical value of the coupling. For N = 100, 



the total planar spin variance reaches its minimum value 
with a coupling value of Ng/ k = —2.034 = N/ (^Jx^, as 
expected. At the minimum variance, we find numerically 
that: 



A^Ji, 



0.1501J= C50. 



(4.4) 



This is in excellent agreement with equation (I2.15P , which 
gives Cgg^ = 0.1499 J. Similar good agreement is ob- 
tained for the calculated values of (^^x'^ , A^J, compared 

with the asymptotic equations p.Sp . (|3.14p . and (|3.15p - 
apart from corrections of order 1/J^. 

In summary, the ground state of a two- well BEC is not 
only a planar squeezed state: at critical coupling it gives 
the exact solution to the minimum variance. However, 
while spin-squeezing has been observed experimentally 
[^-Q , indicating that detection at the quantum shot-noise 
level is technically feasible, existing experiments used a 
BEC with repulsive interactions. To obtain a PQS state 
would require an experiment using attractive interac- 
tions, as found in '^^K, for example [28]. We finally note 
that Einstein-Podolsky-Rosen entanglement and macro- 
scopic superpositions for BEC states have been topics of 
recent interest [22, 25-2% It is intriguing that the sim- 
plest physical route towards generating planar-squeezed 
states also displays other important phy sical properties, 
including macroscopic entanglement |22| . In this case we 
emphasize that the entanglement is found between the 
two underlying boson modes which are used to construct 
the composite spin operators. 



B. Applications to Interferometry 

Due to their interactions with magnetic and grav- 
itational fields, cold atomic sensors have been useful 
for ultra-sensitive magnetometers [29'| and gravimeters 
[30l[3]|. In this section, we analyze the possible applica- 
tions of planar spin squeezing in atom interferometry (3^ . 
In order to improve the performance of this type of sen- 
sor it is important to combine both relatively large atom 
numbers [33] and quantum noise reduction^Q. Consider 
the effect of the spin mapping to a pair of boson modes 
which are subsequently passed through a beam splitter 
(for external degrees of freedom) or microwave rotation 
(in the case of internal degrees of freedom) [1,, Ji] . We use 
the mapping defined in the previous section in equation 
(j4.2p . with a planar quantum squeezed input. 

Next, we consider the effects of an idealized phase- 
measurement. This is shown in schematic form in FigU) 
Two bosonic modes a, b prepared in an entangled PQS 
state, are passed through phase-shifters and a four-port 
50 : 50 beam-splitter. Here <j) is an unknown phase-shift 
to be measured, while 6* is a reference phase-shift. The 
number difference between the output ports, A^"*" — N~ 
is measured, and gives information about the unknown 
phase. 
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Figure 4. (Color online) Schematic diagram of ideal interfer- 
ometric measurements. Phase measurements can be reduced 
in the simplest case to an output measurement of number dif- 
ferences from a beam-splitter, where (j) is the unknown phase- 
shift and S is a known reference phase-shift. 



This has the utility of allowing much lower atom num- 
bers and therefore lower atomic density in atom interfer- 
ometry at a given phase sensitivity. Since the atom den- 
sity is limited by two-body and three-body losses, this is 
a very significant practical advantage. To give an exam- 
ple, a PQS interferometer with 10^ atoms has a phase 
sensitivity of order 10~^ in a single phase measurement. 
To achieve this result with conventional coherent inter- 
ferometry would require 10^ atoms, which implies 100 
times greater density for the same geometry. 



V. PLANAR SQUEEZING ENTANGLEMENT 
CRITERIA 



After the phase shifters, but before the beam splitter, 
the collective spin components in the X — Y plane are: 

Jx{(fy = Jxcos {(f> — 9) + Jysin {(j) — 9) , 
J^((/,) = jycos {(t)-9)- Jxsin {(j) - 9) . (4.5) 
After the beam splitter, we define two output operators 



c± = ae'^ ± be'^ /\/2 



(4.6) 



Calculating the phase-sensitive output number differ- 



ence, N = 



N+ -N- 



and its derivative in terms of 

the measured phase, we find that the measured phase 
uncertainty is: 



A0 = 



v/AW(0) 
\dU/d(i)\ 



(4.7) 



The phase noise in a single measurement in terms of 
the initial spin variances, in a state where (JxJy) — 0, 
is therefore: 



A(f> ^ 



9) + A^JYi 



\Jycos I 



Jxsin ( 



Clearly, in any phase measurement one must avoid insen- 
sitive regimes of the interferometer near the fringe peaks 
where Jy{4>) = dJ\f/d(f> ^ 0. We see that the quantum 
noise in each measurement is bounded above since: 



< A^Jn 



(4.8) 



This shows clearly that the planar variance is an up- 
per bound to the total quantum noise on the interfer- 
ometer number difference output. Squeezing this up- 
per bound below the standard quantum limit is vital 
under conditions where the the measured phase is com- 
pletely unknown. Planar quantum squeezed states there- 
fore can provide useful noise reductions over a large range 
of unknown phases. In the best case, one can achieve 
A^J|| ^ Cj ^ J^/^, which leads to a phase uncertainty 
of: 



A0 oc 



(4.9) 



In earlier works, uncertainty relations have been used 
extensively for the derivation of criteria to detect entan- 
glement ,the EPR paradox and mesoscopic superposi- 
tions ^ H, [H, [13, [M^HI ■ However, neither the Heisen- 
berg [37 1 nor the related Sorenson-Molmer [39l . [4^ in- 
equalities can be used to detect the entanglement of a 
very important subclass of states that have (J) = 0. This 
includes the maximally entangled states which are eigen- 
states of Jz having (Jz) — 0. In these cases, that give 
rise to the classic violations of local realism studied origi- 
nally by Bell [47, 48], the variance of the remaining spins 
are still constrained, because no simultaneous eigenstates 
exist. 

It is this situation our more general two-component un- 
certainty relation involving both variances is extremely 
useful [7]. This two-component criterion can also detect 
the entanglement of the Bell-type maximally entangled 
states, only requiring the measurement of two orthogonal 
spin components [Ti], as well as determining the possible 
existence of multi-particle genuine entanglement. A fur- 
ther application of the planar uncertainty relation equa- 
tion p.ip is that it provides a means of witnessing mul- 
tipartite entanglement between macroscopic spins. 

Hofmann and Takeuchi ff\ have proved that any uncer- 
tainty relation of the type A^A > Ua, where the system 
is labelled A and A is a vector of observables for that sys- 
tem, can be used to define a criterion for entanglement. 
Here the limit U is generically defined as the absolute 
minimum of the uncertainty sum for any quantum state. 
The planar uncertainty equation (jl.ip is of this form. If 
systems A and B are separable, then it is always true 
that 



A^ (A + B) >Ua + Ub. 



(5.1) 



The violation of this uncertainty bound is then a proof 
of entanglement. This result may be generalized to mul- 
tipartite systems consisting of N distinct locations, as 
shown by Toth 0, [1 HO, [ig . Consider N sites of spin J 
particles. Assuming absolute separability, we write the 
total density matrix p as a probabilistic sum of product 
density matrices at site k, occurring with probability 
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R 

Defining collective spin operators as 

N 



jtotal \ ^ jk 



(5.2) 



(5.3) 



fc=i 



where Ck^i = ±1, the expression for the sum of the vari- 
ances is 



N 



A^Jp'>^P«^(A24)«. 



(5.4) 



R k=l 



Using the two-component uncertainty equation (jl.ip . this 
leads to the separability condition for 5*2 , the sum of the 
two-component relative variances: 



2 ftotal 



total > ^ 



(5.5) 



If violated, this will imply an entanglement between some 
of the sites. When three variances are measurable, there 
are similar conditions known based on three- variance un- 
certainties. 

These relations are useful for identifying the entangle- 
ment of maximally entangled states where (J) = 0. For 
example, the Bell singlet-state J — 1/2 for which the 
spins are anti-correlated {^J*°*°-^ = 0) gives the total un- 
certainty as zero (i.e., 5*2 = 0) provided = = = 1. 
The many body singlet state in which pairs of parti- 
cles are in singlet states gives zero total uncertainties for 
larger A'' [i^ . Here we consider the maximally entangled 
states for fixed J of form (here | J, m) are the eigenstates 
of Jz) 



I J, TV) 



1 



M 



V2J+1 



(5.6) 



m— — J 



where spins are correlated. For N — 2, the sum of the 

two variances A ( — J^) + A (^Jy + Jy) is always 
zero, for any J. Singlet states of total spin zero for d- 
level systems have been presented in |50]. We denote the 
singlet state (which has zero total spin) obtained from N 
particles of spin J as | J, N)s- 

The two variance sum ^2 is zero for the anti-correlated 
singlet state. To examine the sensitivity of the two vari- 
ance criterion to noise, we reduce entanglement by con- 
sidering the mixed state of the type considered by Werner 

P = pJj + (1 - Pn) \J, N)ss{J, N\ , (5.7) 
where p„ gives the relative contribution of the white noise 



term represented by Ij 



1/ (2 J + 1) 



which is a ro- 



tationally symmetric, uncorrelated state proportional to 




Figure 5. (Color online) Plot of the multipartite entanglement 
criteria. The sloping, bold lines give the normalized sum of 
the two spin variances {S2 /NJ) for the A'^ particle singlet state 
\J,N)s mixed with the maximally noisy state, as a function 
of noise probability p„. Entanglement is confirmed when the 
bold line for each J falls below the corresponding horizontal 
line giving Cj/J. From the bottom, the lines correspond to 
J = 1/2,1,2,10. 



the identity matrix at each site. Since the singlet state is 
perfectly correlated, and the Werner state is completely 
isotropic, one can show that the uncertainties for the mul- 
tipartite case are entirely due to the white noise terms. 
From equation (|5.5p . which requires only two measure- 
ment settings, the condition for detecting entanglement 
is given by: 



S2^^J{J+l)Pn<NCj. 



(5.8) 



which gives the bound of noise for detecting entanglement 
of: 



Pn < 



3C,7 



2J(J-|- 1) 



(5.9) 



If all three spin measurements are available, a three spin 
measurement is even more sensitive to entanglement for 
this isotropic case. However, in many cases all three com- 
ponents are not measurable, or may not all have strong 
correlations. 

Figure [5] illustrates the purity required in order 
to demonstrate entanglement between spin J particles 
(sites) using the two variance equation (|5.5p . Entangle- 
ment is detected when the normalized two variance sum 
5*2 is below the line indicating the value of Cj/J, i.e., 
when S2/NJ < Cj/J. 



VI. OUTLOOK 

In conclusion, we have derived an uncertainty relation 
for the planar sum of the variances in two orthogonal 
spin directions for systems of fixed total spin. The lower 
bound varies asymptotically as J^^"^. We have shown that 
this planar local uncertainty relation can be readily sat- 
urated with macroscopic planar squeezed states at large 
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spin. In addition, a practical technique for generating 
these states is proposed, employing the ground state of a 
two- well Bose condensate with attractive interactions. 

Such planar squeezed states have the feature that they 
are able to minimize phase measurement noise over a 
wide range of unknown phase angles, in an interferomet- 
ric measurement. Since these states are readily obtain- 
able as ground states of physically relevant Hamiltonians 
in a two-mode Bose-Einstein condensate, it appears fea- 
sible to generate and demonstrate these features in labo- 
ratory measurements, either using ground state prepara- 
tion in an attractive BEC [22| . or using dynamical tech- 
niques with repulsive interactions. 

Criteria for detecting entanglement between multiple 
spin J systems using only two component measurements 
can be derived from this. The two spin component entan- 



glement criterion is likely to have important applications 
where noise or measurement is asymmetric, so that all 
three components cannot be measured. There are other 
generalizations possible if one constrains {Jz) to have a 
finite value, giving an inequality involving Cj {Jz) ; these 
will be treated elsewhere. We expect similar relations to 
occur for other continuous groups. 
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